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APPLICATION  OF  THE  OPERATOR  EXPANSION  METHOD  TO 
REALISTIC  TWO-DIMENSIONAL  (2D)  SEA  SURFACES 

JOHN  DUDDHRLHY  and  RICHARD  S.  KEIPHiR 
Naval  Research  I -ahoratory 
Siennis  Space  Ccnicr,  MS.  USA  39529-5004 


SUMMARY 

Tlie  ofKrator  expansion  nicihod  has  been  shown  lo  model  quickly  and  accurately  acoustic 
scattering  from  surfaces  with  roughness  similar  to  that  of  the  oceans.  In  this  paper,  the 
monostatic  backscattcr  is  computed  using  the  operator  expansion  method  over  an  casemble  of 
modeled  2D  deterministic  st'a  surfaces.  In  particular,  low  gra/Jng  angle  azimuthal  dependence 
of  tlic  backscattcr  is  explored. 


1.  INTRODUCTION 

Tlie  operator  expansion  method  for  quickly  .solving  die  Helmholtz  Integral  Equation  for  the 
scattered  field  was  outlined  by  Milder •.  Tliis  method  finds  the  scattered  fields  at  an  interface  by  a 
Taylor  series  approximation  with  respect  to  I’  c  surface  height  based  on  the  filtered  FFT  of  the 
incident  fields.  Since  the  operator  expansion  mt  /w'd  is  driven  by  FFTs  rather  than  by  either  direct  or 
indirect  matrix  inversion,  it  is  much  quicker  tlian  exact  integral  equation  methods.  Furlhcnnorc  as  the 
number  of  field  points  increases,  the  calculation  time  increases  at  a  much  slower  rate  for  the  operator 
expansion  meth^  than  with  the  exact  integral  cqua  ion  methods.  The  operator  expansion  method  has 
also  been  shown  to  be  sufficiently  accurate  for  regime  of  interest  here.  Kaezkowski  and  Thorsos 
demonstrated  that  this  method  spans  the  regimes  of  validity  of  both  the  Kirchhoff  and  small 
perturbation  approximation^.  Furtlicr  confirmation  by  the  authors^  has  shown  that  third  order  and 
higher  terms  of  the  symmetric  operator  expansion  arc  unnecessary  for  modeling  rough  sea  surfaces 
reflection  for  the  source  frequency  range  considered.  Here  we  will  show  the  average  monostatic 
backscattcring  from  rough  surfaces  representative  of  a  surface  wave  spectrum  encountered  in  the  Gulf 
of  Alaska  wiiliout  any  consideration  of  po.s.sible  bubble  fields.  Our  goal  in  accomplishing  this  is  three 
fold.  1)  We  wish  to  show  that  a  wide  range  of  scattering  problem  that  were  previously  too  large  to 
compute  can  now  be  .solved.  2)  We  hope  to  gain  some  sort  of  insight  into  the  scattering  behavior  of 
two  dimensional  random  surface  fields.  3)  By  examining  simulated  sea  surface  scattering  without  a 
consideration  of  bubble  fields,  we  wish  to  learn  more  about  the  contributions  to  observed  sea 
scattering  made  by  both  surface  scatiing  and  bubble  scattering. 

2.  SPECTRUM  AND  SURFACE  GENERATION 

If  a  survey  of  the  literature  were  conducted  on  the  subject  of  modeling  the  acoustic  scattering 
from  the  sea  surface  roughness,  it  would  reveal  that  almost  all  of  the  analysis  to  date  involves  surface 
spectra  due  to  tcmporaliy-stationary,  spatially-uniform  forcing  conditions.  If  two-dimensional  (2D) 
sea  surfaces  are  considered,  a  simple  azimuthal  dependence,  perhaps  not  even  wavenumber 
dependent,  is  typically  assumed.  In  reality,  the  .sea  surface  is  influenced  by  nearby  as  well  as  distant 
forcing  phenomena,  evolving  over  time  in  a  complex  manner.  In  order  to  examine  acoustic  scattering 
from  sea  surfaces  that  arc  due  to  temporally  non-.staiionary,  spatially  non-uniform  generating  winds, 
i.c.,  more  realistic  sea  surfaces,  the  DWAVE  Model  has  been  employed.  This  modcl'^  numerically 
calculates  the  directional  wave  spectrum  at  a  given  geographical  location  ba.scd  on  the  time-varying 
wind  field  and  other  environmental  information  spccifial  over  an  extended  region.  Tliis  2D  numerical 
spectrum  is  then  used  in  tlie  filtering  process  that  is  described  below  in  order  to  generate  cascmbles  of 
scattering  surfaces. 

The  particular  spectrum  considered  in  this  paper  occurred  during  the  spring  of  1990  in  the 
Gulf  of  Alaska.  During  the  hours  prior  lo  this  event,  winds  in  the  area  ranged  from  8-15  m/s.  The 
rms.  roughness  of  the  .seas  was  .296  meters  at  this  time.  A  contour  plot  of  the  surface  roughness 
spectrum  is  shown  in  figure  1.  In  order  lo  generate  realizations  having  this  spectrum,  a  2D  field  of 
Gaussian  deviates,  were  Fourier  transformed,  filtered  in  the  transform  domain,  and  Inverse 
transformed.  In  this  way,  the  2D  array  of  random  surface  heights  arc  correlated  to  have  the  desired 
roughne.ss  spectrum.  For  the  cases  considered  herein,  the  surface  spectrum  was  sampled  512  times  in 

both  the  kx  and  ky  at  a  rate  of  Akx  =  Aky  =  .(H  r:id/m.  1his  resulted  in  malizations  defined  at  intervals 
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Figure  1  Logarilhmic  contour  plot  of  the  chosen  sea  surface  roughness  spectrum. 
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Figure  2.  Backscaller  intensity  plot  versus  azimuth  Figure  3.  50  by  50  point  sample  of  the  first 
for  tlic  surface  shown  In  figure  2.  .surface  realization  of  the  cnsomhlc. 


\x  -  Ay  -  .MY!  iiiul  covering  an  area  ol  neatly  25,0(X)  stiuare  meters.  For  Hie  calculation  of-the 
staiisiics  ol  (he  .scallered  field  .'i  '.sea  surface  rcali/alions  were  used. 
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For  this  data  set  we  examined  the  a/imuihal(((>)  dependence  of  tlic  inonostatic  backscatter. 

I  lolding  Hie  incident  grazing  angle  O'f^'  fixed  at  ten  degrees  backscatter  was  calculated  for  every  two 
degrees  in  azimuili  llie  2fi()  tF/.  incident  pre.ssure  wave  fields  calculated  at  each  sample  point  for  each 
incident  plii  were: 


pine(p)  _  expl-iVr  - 


4((x^  +  y2  _  f  coi(0‘*'‘-')) 

_ 
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Wliere  1 .  =  the  surface  sample's  extent  in  ilie  x  and  y  directions,  the  wavelength  is  6  meters,  and  the 
origin  is  set  at  ilic  center  of  tlie  surface  sample.  Dack.scaiier  in  dB  is  calculated  relative  to  tlie  far  field 
scaitering  intensity  given  by  a  scattercr  that  is  non  directional  in  plii  plane  of  incidence.  Tliis  can  be 
expressed  fornuilaically  in  tlie  following  way: 
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Wliere  is  the  far  field  normalized  pressure  field  in  the  given  direction,  P*  is  the  pressure  at  each 
of  the  n  surface  points,  and  m  is  (lie  numlicr  of  ihela  samples  taken  in  the  phi  plane.  Figure  2  shows 
the  back.scaiter  of  tlie  above  acoustic  beams  from  a  single  surface  realization  which  has  been  partially 
illiistraied  in  Figure  3.  Figure  4  shows  the  average  backscatter  in  phi  for  the  50  surface  ensemble  as 
well  as  tlie  90%  confidence  range  in  dD.  Calculated  but  not  shown  is  the  perturbation  theory  results 
for  tlie  same  spectrum.  Perturbation  theory,  though  lacking  details,  gives  the  same  general  shape  in 
azimuth.  In  particular  the  null  around  90  degrees  in  azimuth  in  figure  6  is  in  agreement  with  first  order 
penurbaiion  theory.  All  tlie  above  calculatioas  were  performed  on  a  Cray  Research  Corporation  Y-MP 
eiglil  vector  processing  supercomputer. 


Phi  In  Degrees 

I  igiiie  4.  Average  Baclcscatler  of  50  clement  eitseiiihle  ol  surface  realiz.atioas  of  the  spectrum  in 
liguie  2.  Backscatter  given  in  dlJ  relative  to  isotropic  scattering. 
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4.  CONCLUSIONS 

llic  s  'iiiiioiric  opcralor  expansion  enables  one  to  iiiiickly  calculate  the  acoustical  energy 
•jrlace  seaiicfiiig  with  less  than  one  pan  in  a  thousand  error  for  surfaces  ttxi  large  to  compute  using 
olhci  sir.iilarly  accurate  metJiods.  For  example  a  single  backscatter  geometry  takes  approximately  45 
seconds  ttt  solve  using  the  operator  expansion  mcUiod  versus  over  a  hundred  days  for  a  similarly 
accurate  Iterative  Moment  Method  solution  for  the  same  geometry.  Furthermore,  while  the 
Periurbatittn  Tlicory  method  produces  a  quick  look  at  the  approximate  ensemble  average  backscatter 
from  a  given  surface  wave  spectrum,  the  symmetric  operator  expansion  method  has  the  advantages 
over  Penurbation  Theory  that  all  scattering  geometries  in  the  half  space  can  be  calculated,  non-plane 
wave  incident  fields  can  be  introduced,  more  statistical  data  can  be  produced,  and  grater  accuracy  can 
be  achieved. 

These  calculations  show  a  strong  dependency  of  the  backscatter  intensity  on  the  azimuthal 
angle  for  all  surfaces  generated  in  the  data  set..  This  dependency  is  also  visible  in  all  other  calculations 
performed  on  modeled  Gulf  of  Alaska  surface  spectra.  This  observation  implies  that  the  rough  surface 
backscatter  status  as  cither  possibly  contributory  or  non-contributory  to  measured  ocean  baclcscatter 
depending  on  whctlter  a  measured  phi  angular  backscatter  dependence  is  respectively  found  or  not 
found.  Unfortunately,  there  is  no  published  data  found  by  the  authors  tltat  conclusively  answers  this  _ 
question.  Hopefully  experiments  will  be  performed  or  data  will  be  published  that  could  conclusively 
answer  this  question. 
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